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1. Introduction 

Universality and scaling are two important concepts in the theory of critical phenomena 
|I|, [J and the Ising model P has been a model widely used in such studies. Recently, 
exact universal amplitude ratios and finite-size corrections to scaling in critical Ising 
model on planar lattices have received much attention [0 |5], 0, H |H], HH HI 



This may be due to the fact that the hypothesis of universality leads naturally to the 
consideration of universal critical amplitudes and amplitude combinations |13| , and for 
the comparison between experiment and theory in relation to scaling and universality, 
it is often a more rigorous test to use amplitude relations rather than critical exponent 
values. Moreover, it is also well known that the finite-size scaling functions depend 



on the boundary conditions [14|, and there has been considerable recent interest in 



studying lattice model with various boundary conditions |T5], [16], [IT], [18], [19], |20|, [21 
The study of exact universal amplitude ratios and finite-size corrections to scaling in 
critical Ising model is usually based on the analytical solutions of the model on finite 
lattices. Although the exact solution of the Ising model on M x N square (sq) lattice 
had been obtained long time ago |22| , and the exact expression of the partition function 



of the Ising model on M x N plane triangular (pt) lattice has been obtained by lattice 
field theories recently |[23[| , there is still no published results for the exact solutions of the 



Ising model on M x N pt and honeycomb (he) lattices with periodic-aperiodic boundary 
conditions. The purpose of this paper is to fill this gap. In the present paper we use 
the Grassmann path integral to calculate exact partition functions of the Ising model 
on M x N sq, pt and he lattices with periodic-periodic (pp), periodic-antiperiodic (pa), 
ant iperio die-periodic (ap) and antiperiodic-antiperiodic (aa) boundary conditions. The 
partition functions are used to calculate and plot the specific heat, C/ks, as a function 
of the temperature, 9 = ksT/J. We find that for the N x N sq lattice, Cjk^ for pa 
and ap boundary conditions are different from those for aa boundary conditions, but 
for the N x N pt and he lattices, C jks for ap, pa, and aa boundary conditions have the 
same values. Our exact partition functions might also be useful for understanding the 
effects of lattice structures and boundary conditions on critical finite-size corrections of 
the Ising model. 

Two-dimensional Ising model on the sq lattice at vanishing magnetic field was first 
solved by Onsager by the use of Lie algebra [[|. The exact solution he obtained was 
Ising model on an infinite lattice. The original method was rather complicated, and it 
was later improved by Kaufman |22|] who obtained the exact solution of the Ising model 
on a finite torus by using the theory of spinor representation. The successful treatments 
of the two-dimensional Ising model brought the studies of phase transition into the 
modern era. Onsager's solution in one hand showed the previous classical theories 
were unreliable in their quantitative predictions, and on the other hand provided a 
great stimulus to explore the true behaviour near the critical point. After Onsager's 
original solution, many quite different mathematical approaches were developed, but the 
approaches were still complicated. Among them, Schultz, Mattis and Lieb gave explicitly 
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the fermionic treatment in the framework of transfer- matrix formalism fEM], and Kac 



and Ward developed the combinatorial method |25], |26| . Both methods reformulated the 
two-dimensional Ising model as a free-fermionic field theory in terms of anticommuting 
Grassmann variables, which enclosed the fact that the Ising model on two dimensional 
regular lattices may be viewed as free-fermionic theory The other alternative method 
in literature was the Pfaffian representation, which was introduced by Kasteleyn |27| to 
translate Ising spins into dimers that can be reduced to some Pfaffian ||28|| . Stephenson 
has used the Pfaffian representation to solve the Ising model on the pt lattice, but 
the solution was restricted to 6L x 6L lattice due to its 6 x 6 basic nonvanishing 
matrix elements and was exact only in the limit of L —>■ oo [p9fl . Recently, by using 
the connections between Pfaffian, dimer and Ising model, Nash and O'Connor have 
obtained the exact expression of the partition function of the pt lattice Ising model on 
a finite torus [23||. They first employed the lattice field theories to obtain the exact 
partition function of the Gaussian model, and then established the exact expression of 
the partition function of the pt lattice Ising model from the analysis of the appropriate 
lattice determinants and the parameterization according to the results in p9| . 



On the other hand, in view of the simplifying the approach, a remarkable progress 
was achieved by Plechko who modified the traditional fermionic interpretation and 
introduced a nonstandard approach By the use of this approach, Plechko himself 
has not only rederived Onsager's and Kaufman's results in a relatively simple way 
but also obtained the partition functions of a class of triangular type decorated 



lattices pl| , and a triangular lattice net with holes f3"2| . Quite recently, by using the 
same approach, Wu at al. have obtained the M x iV sq lattice Ising model with 
periodic-aperiodic boundary condition ||, and Liaw at al. have successfully solved 
triangular and hexagonal lattices on a cylinder geometry (M x oo) with periodic and 
antiperiodic boundary condition [33]. This approach is based on the integration over 
the anticommuting Grassmann variables and the mirror-ordered factorization principle 
in two-dimensional density matrix |30, |3l], [32], [33|, and does not involve the traditional 
transfer-matrix or combinatorial considerations. The whole scheme of the method can 
be illustrated schematically as shown below [jBOfl : 

Z = Sp {Z (a)} -> Sp {Z (a \ X )} -> Sp {Z &)} = Z, 

{*) Hx) (x) 

where w Sp" stands for the average over spin variables (a) or Grassmann variables (x)- 
The original partition function Z is expressed purely by spin variables (a) at each 
lattice site. With a set of anticommuting Grassmann variables (x) being introduced to 
factorize the local bond Boltzmann weight such that spin variables are decoupled, the 
partition function passes to a mixed Z (a \ x) representation. Then, by eliminating the 
spin variables in the mixed Z (a\x) representation, the fermionic interpretation Z (x) of 
the two-dimensional Ising model can be obtained, and after carrying out the Grassmann 
integral, the analytic solution for the partition function and free energy can be achieved 
53f. 



31, 32 



In the present paper, we work in this framework to obtain exact partition functions 
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of M x N pt and he lattices with different boundary conditions, including pp, pa, ap and 
aa boundary conditions. We used these results to calculate and plot the specific heat, 
C/ks, as a function of the temperature, = ksT / J . Our results show that for the sq 
lattice, C/ks for pa and ap boundary conditions are different from those for aa boundary 
conditions, but for the pt and he lattices, C/ks for ap, pa, and aa boundary conditions 
have the same values. Beside these analyses, our exact partition functions may also 
be used for understanding the effects of lattice structures and boundary conditions on 
critical properties and critical finite-size corrections of the Ising model. 

This paper is organized as follows. In section 2, we set up a general form of the 
partition function for pt and he lattices. Then, three pairs of conjugate Grassmann 
variables are introduced for a lattice site to factorize the Boltzmann weights, and the 
principle of mirror ordering are used to rearrange the Grassmann factors so we can 
perform the summation over Ising spins to obtain a pure fermionic expression of the 
partition function. In section 3, using the Fourier transform technique we complete 
the integrations over the Grassmann variables to obtain the exact solution of the 
partition function. Then, the solution is subjected to periodic-aperiodic boundary 
conditions, including pp, pa, ap and aa boundary conditions. We further consider 
the shift behaviours of the maximum of the specific heats of these systems in section 4. 
Finally, we discuss some problems for further studies in section 5. 



2. The Partition Function 



Consider Ising ferro magnets on M x N pt and he lattices as shown in figure 1, in which 
the former is considered as a sq lattice with a single second-neighbor interaction, and 
the latter contains an inner spin in each lattice cell. The corresponding Hamiltonians, 
respectively, read as 



M N 



Ht — — (JlO~ mn O~ m +ln + ^^mn^ran+l + J m+lnO~ mn+l) , 



(1) 



and 



rn=l n=l 



M N 



H h = 

m=l n=l 

where Jj with % = 1,2,3 are the coupling constants (Jj > for ferromagnetic lattices), 
a mn = ±1 is the Ising spin located at the site (m, n), and <7o denotes the inner Ising 
spin in he lattice. Using the identity of the Boltzmann weight, 

exp (PJi(T^a u ) = cosh (/3Jj) [1 + tanh (/3Jj) a^} , (3) 

j3 = (ksT) -1 , and performing the sum over <t , the partition functions of two lattices 
can be formulated in a single three spin-polynomial representation, 



n b 



jQcosh (/3 Ji] 



.1=1 
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{M N ~\ 
n n ( a ° + 
m=l n=l J 

where iV, is the numbers of lattice sites (N s = MN for sq and pt lattice, N s = 2MN for 
he lattice) and n& is the number of bonds per lattice cell (n& = 2 for sq lattice, n& = 3 
for pt and he lattice), symbol vv Sp " stands for spin average defined by 

(a) 

SpH = £ £ [•■•]' S p[l] = l, SpN = (5) 

fa) Z (a l= ±l) (<Tl) 

and ctj's are defined as 

«o = 1 + tit 2 t 3 , = ti + t 2 t 3 , = t 2 + hh, = t 3 + tit 2 , (6) 
ti = tanh with i = 1, 2, 3, for pt lattice, and 

= 1, af = tit 3 , of = ht 2 , af = t 2 <3, (7) 

for he lattice. 

To factorize the partition, we rewrite the partition function as 

n b 

Jlcosh {13 Ji 



Z H = 2* 



N s 



_i=l 

r m n 

x Sp < J| J| r (1 + r l°~mn°~m+ln) (1 + r 2°'mn 'mn+l) (1 + r 3 cr m+ln cr 
fa 



mn+1) ( i 

m=l n=l 



where r» with i = 0, 1,2,3 vary from one lattice to the other, and are related to ojj's 
from 

a = r (1 + nr 2 r 3 ) , a x = r (n + r 2 r 3 ) , a 2 = r (r 2 + nr 3 ) ,03 = ^0 (r 3 + rir 2 ) . (9) 

For pt lattice, the relation between rj and ij is trivial, i.e. ro = 1 and r« = £j, but for he 
lattice, the relation is nontrivial and is determined by equations (|7|) and (Q). 

It is more convenient to define the generalized reduced partition function as 

Q = r* IN Q, (10) 

with 

M N 

Q=\\ n s P r l°~mnO~m+ln) (1 + r 2® mnO~ mn+\) (1 + r 3°~mn+lO~m+ln)] ■ (H) 

TO =1 n=l fa«») 

To construct the fermionic representation of the generalized partition function, 
we associate each lattice site (m, n) with three pairs of conjugate Grassmann variables, 
{a mn , a* mn ; b mn , b* mn ; c mn , c* mn } G X- All of these Grassmann variables are anticommuting, 



and their square are zero. Their integral obeys the basic rules |34| 



dx = 0, J d X -X=h (12) 

d x -{i(x + v)= I d x -n( x ), (13) 
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for an arbitrary vector rj with anticommuting components, and there is the relation 

1 + na^ = J d X *d X e xx * (1 + (1 + nx*^) • (14) 



Using these Grassmann variables, we can rewrite the reduced partition function aspO 

M N 

Q ] [ \\ Sp Sp (A mn A^ +ln B mn B^ n+1 C mn+ iC^ +ln ) 
where vx Sp " stands for the averaging with Gaussian weight 



(15) 



Sp [•■■]= / d X *d X ie™* [■■■}, (16) 

with the rules 

Sp [XiX*] = - Sp [XiXi] = 1, (17) 
(Xi) (x») 

Sp [x»] = Sp [*?] = 0, (18) 

(Xi) (x0 

and the Grassmann factors, A, A*, B, B*, C, and C*, are defined as 

^mri 1 "I - a mn o~ ran ^ A mn 1 -|- r\a m _^ n o~ mn} (19) 

B m n = 1 4~ b rnn O' mn , B mn = 1 + T ^ mn _^(J mn \ (20) 
Cmn 1 "I - *--mn— l""mn> ^ mn 1 '~3'-' JTl _i n O~rfin- (21) 

In this way, a Boltzmann weight is decoupled to the product of two factors of separated 
spins. 

For simplicity, we express the reduced partition function as 

{M N ~\ 
n n *^*2»*s» » (22) 
m=l n=l ) 

where #^ n , and are defined by 

^mn = Sp (A mn y4 m+ln ) , (23) 

(°~m?i ) 

^mn = S P { B rnnB* mn+ ^) , (24) 

(°~m?i ) 

^mn = Sp (C mn+ iC^ +lrl ) . (25) 

We first treat the boundary weight and consider periodic boundary condition in 
both directions: 

^Mn = Sp [(1 + a Mn a Mn ) (1 + ria^ n cr M +in)] 

(o"Mn) 

= Sp [(1 + rial n ai n ) (1 + a tfn ffMn)] 

(o"Mn) 

= Sp K n AMn), (26) 

(o"Mn) 
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which implies 



A 0n 



Similarly, from 



we have 



* mN 



* Mn 



^ mN 



h* 



On. 



h Mn- 



Sp (B mN B* mN+] ) — Sp (B* ml B mN ) , 

( CT m]v) (fjVn) 
Sp (CA/n+lCIf+ln.) = Sp (Cl n C Mn+1 ) , 

("Mn) {"Mn) 

Sp (CmN+lCm+m) = Sp (C»nlCm+lAf) 

( CT mJv) ( CT mJv) 



u mNi 



fC 



(27) 

(28) 
(29) 
(30) 

(31) 
(32) 
(33) 



Since c m0 = c m N, need not to be treated as a boundary weight, and only ^ Mn , 

^ mN and ^ Mn should be considered. However, this situation becomes ambiguous when 
we take Fourier transform of these Grassmann variables with single set of exponential 

). Because the Fourier exponential factors are associated 



factors in equations (pq) and 
with directions in M and N, the sign factor in front of h* mN takes effects simultaneously 
on b* mN and c m ^. Therefore, the real situation is that instead of the relation in equation 
(j33|), we must take 

CmO = —CmN- (34) 

A self-consistent way to assign a minus sign to c m o and obtain the relation in equation 
([34]) is interchanging C m ± in equation ( p0|) with another Grassmann factor. An 
equivalent but more convenient approach is to consider the rearrangement of B* mX in 
the boundary weight together with the rearrangement of C m \ in the reduced partition 
function. To see this, we express the reduced partition function as 



'M-l N 



TT TT v A v c 

1 JL JL A mn 1 



Q = Sp < Sp 

(a,b,c) [ (<x) Vm=1 n=1 

with the boundary weight ^> b 



Sp 

(a,b,c) 



B ■ 



M N-l 
vm=l n=l 



R* 
mn D ran+l 



(35) 




Ml 



M 

n 

m=l 



(36) 



M-l N 



m=l n=l 



n n *i»*s« = n n a ^°- 

m=l n=l 
M-l / N 

= Y\_ I IT CmnAmn J C-ml I \\ ^m+ln^m+ln 



(37) 



m=l 



vn=2 



vn=l 
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Here, arrows have been used to indicate the orders of the products in m and n. When 



A' 



N 



we move B ml from left of Yl C* n A\ n to right of Yl Cin^in m equation 

n=l n=l 

N n 



)> B *ml 



passes 2N Grassmann factors, but for moving C m i from right of ]~[ C mn A mn to left of 

n=2 

A m i in equations (|36| ) and (|37|), C m i passes only 2iV — 1 Grassmann factors. Then by 
moving B ml from left to right, and simultaneously moving C m i from right to left, we 
can assign to the Grassmann variable in C m i an additional minus sign comparing with 
the Grassmann variable in B ml , and hence obtain the relation of equation ([34]). 

Accordingly, we interchange B* and C* A* in equation ( pq) to obtain the 
arrangement of C*A*B* according to the identity 



B + (CA) + = - [{CA) + B + + (CA) + B- + (CA)~ B + - {CA)~ B~] 



(3f 



with superscripts + and — being the sign factors in boundary Grassmann factors 



N 



A\ n , B ml and and simultaneously move C m i from right of Yl C mn A mn to left 

n=2 

of A m i in equations ( j3~6|) and (j3~7f). Here we note that the superscripts + and — 
respectively correspond to periodic and antiperiodic boundary condition imposed on the 
spin variables and in turn on the Grassmann variables. Hence, the reduced partition 
function becomes 

1 



Q 



with 



Qi = Sp <^ Sp 

(a,b,c) I (it) 



'M-l 



r 2 



+ Q, 



r 3 



r 4 



M N-l 



n e - e - + i • ^ • 1 1 n 



mn+l 



. m=l 



Vin=l n=l 



and 



7 



Sp < Sp 

(a,b,c) (a) 



M 



M 



©i II b.u ][ 



->mjV 



vm=l 



\m=l 



where we have defined 

N 



N 



P) — TT C A * and ft* — TT 4*~ 



(39) 



(40) 



(41) 



(42) 



71=1 



71=1 



and the boundary conditions r 1; T 2 , T 3 , T 4 are defined as 



Ti = 


( fl 0n — 


* 

~ a Mni 




— h* 
u mNi 


C 0n — 


_C Mn) 5 


(43) 


r 2 = 


l a 0n — 


* 

~ a Mni 






C 0n = 


_C Mll) > 


(44) 


r 3 = 


( a 0n = 


+ a M,n> 


h* — 


— h* 

u mNi 


r* — 


+ C M,n) 5 


(45) 


r 4 = 


( a 0n = 


+ a M,n> 


h* — 


^mNi 


C 0n = 


+ C M,n) • 


(46) 



In this way, the configurations of the reduced partition function can be further 
rearrangement and expressed as 



'M-l 



M 



M 



M N-l 



sp s P h ii e^r; ei 9m mt nn^ 



(a,b,c) (a) 



R* 



x m=l 



\m=l 



\m=l 



m=l n=l 
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{/ M m \ / M n \ / M N-l \ ^ 

n Q*mB* ml Q m n ^ n n B ^ mn+l . (4?) 
\m=l / \m=l / \m=l n=l / J 

To have a complete mirror-ordered form, we have to rearrange the terms in the last 
two brackets. To achieve this, first we note that 

> 

M /N-l n \ / M m \ M N-l 



<5 7 = ^ Sp ^ Sp ^ Y\ Qm^ml ( Y\. CmnAmn ) C mN A mN ( £? mM J B mn B, 



mn+1 

m=l \n=l / \m=l / m=l n=l 



M /N-l n \ > / M m_ 

Sp Sp < ] [ Q* m I J J B mn C mn A mn B mn J B mN C m NA mN I J [ B mN J ^ . (48) 

(«A«0 (<r) | m=1 \ n=1 / \m=l 

The boundary term of with m = M, denoted by T, can be formulated 

/N-l n \ 

T = @M I | [ B* Mn C MnA\JnB Mn J B* MN C M N Am N B \J N 



,n=l 



N 



W C Mn A* Mn I I Y\ B* Mn C Mn A Mn B Mn 

\n=l J \n=l 

N 

= \\C Mn A* M Ji Mn C Mn A Mn B Mn , (49) 

n=l 

due to the fact that Sp [C Mn A Mn B Mn CMnAM n BMn] for a given n is a commutable 

(&mn ) 

object. By continuing such construction from m = M down to m — 1, we can obtain 
the expression 

{M N \ 
W \\ Sp [C mn A mn B* mn C mn A mn B mn ] > . (50) 
m =l n =l J 

For this partition function, the factors containing the same spin are grouped 
together and we can perform the average over spins. As a result, we have 

/M N / M N \ 

II II da *mn da m ndb* mn db mn dc* mn dc mn exp ( F mn J , (51) 

m=l ra=l \m=l n=l / 



with 



Fmn — a mn a ' mn + b mn b mn + C mn C mn 
+ r l r 3 C m-ln a m-ln 



+ (^C_i n + ria^_ ln ) r 2 b* mn _ x 

+ ( r 3 c m -l n + r l a m-ln r ^mn-l) c mn-l 

+ ( r 3 c m -ln "I" r 'l a m-ln T '^mn-1 c mn-l) Oran 

+ { r z c m -in + r i a * m -in + r 2b* mn -! + c m „_i + a mn ) b mn . (52) 
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Since there is no mix on a mn and b mn , the integral in the above expression can be 
simplified by integrating out the a mn and b mn fields by means of the identity 

J dbda exp {Xab + aL + L'b) = A exp (X^LL 1 ) , (53) 

where a, b are Grassmann variables, L, V are linear fermionic forms independent of a, 
b and A is a parameter. The result then becomes 

M N / M N 



Q-y /II J^J dg mn dg mn dc mn dc mn 

m=l n=l \m=l n=l 



(54) 



with 

^mn — c mnC mn T~ 9mn9mn 

+ ?" irzc* m _ ln g m -i n 

m—ln '"lfl'm— Ira) ^'29mn—\ 
+ ( r 3 C m-ln + r l9m-ln — r 29 mn -l) C mn -1 

~~ ( r 3 C m-lra + r l9m-ln ~ r 29mn~l + c mra-l) (?mn + 3 ran ) j (55) 

where we have changed the notations for the fields by (a^ n ,6^ n ) — > (<? mn , ~9 mn )- This 
is the pure fermionic representation of the reduced partition function. 

3. Exact Solution 

Next, to carry out the integration, we have to use the technique of Fourier transform to 
treat the Grassmann variables which mix together with the variables at different sites. 
The Fourier transformation is defined as 

1 M-l N-l 

X ™ = E E X pq e-^e-^, (56) 

V p=0 q=0 

and 

M-l N-l 

X* = = Wl* e ^ mp e^ n9 (57) 

where the variables X mn and X^ n denotes one of the variables {c mn , g mn } and {c* mn1 g mn } 
respectively. 

After performing the Fourier transformation, the partition function becomes 

Af-lJV-l „ 

Q 7 =IIII dV pq exp(H pq ), (58) 

p=0 q=0 J 

with the measure dV pq defined as 

dV pq = dg* q dg pq dc* pq dc pq , (59) 
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H 



,„. - y l -r 3 e l M p e l N q ) c pq c pq 



+ [r 2 - e i2 * q I c vq g; + r 3 I n - e I 



(l + ri e c pq g M - P N- g - r 3 e ^ (I + r 2 e^ q ) c* pq g* M _ pN _ 

r 2e -*N-i - r x r 2 e i ^ lV e~ i ^ q ) g pq g* 



+ 1 — rie 



1 AI P 



; 2ir , 



Jpq 



- r 1 e t MPg pg g M ^ pN _ q + r 2 e " L ^g; q gh- P N- q 



(60) 



Because H pq contains not only the variables, X pq and X* q , but also the variables 
X 



M-pN-q 



and X 



M-pN-qi 



instead of calculating Q-y it is easier to calculate Q given 



by 



Here H 



M-pN-q 



M-17V-1 „ 

<3?=nn/ dV pq dV M - P N- q exp (H pq + H* M „ pN _ q ) . (61) 

p=0 g=0 ^ 

can be obtained from H pq by replacing p by M — p and q by iV — q for 
the Grassmann variables and replacing the coefficient in front of Grassmann variables 
by its complex conjugate. Completing the integration yields 

Q 7 = n if U - M cos ^ - A 2 cos ^ - A 3 cos - 1 , (62) 



p=0 g=0 



with 



An 



M 



2 2 2 2 



iV 



Ai = 2 (a ai - a 2 a 3 ) , 
A 2 = 2 (a a 2 - «ia 3 ) , 
A 3 = 2 (a a 3 - a x a 2 ) , 

where a , ai, a 2 , and a 3 are given by equations @ and (0) for pt and he lattices, 
respectively. 



(63) 
(64) 
(65) 
(66) 



3.1. Periodic-periodic boundary condition 



According to equation (|39"D, the reduced partition function for ferromagnetic lattices 
with pp boundary condition is 



Q 



pp 



2' 2 2> u u, 2 



sgn 



6-9 c 

~o7 



no 



(67) 



1/2 



.(68) 



where the superscript p refers to periodic boundary condition and 

M-l N-l 

^=U II [Ao-Acos^-A 2 cos^-A 3 cos(^-^) 

p=0 g=0 

The sign factor in front of the last term is a result of the standard consideration of the 
Grassmann integral over the zero-mode variable p = q = for ferromagnetic couplings 
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35| . When the integral of equation (^) is carried out, it is always positive, but 
this is not the case for equation (|58|). There are unpaired terms from zero-mode in 
equation ( p8|) under various boundary conditions and they contribute a sign factor to 
Q4 for < ij < 1. The partition function for pp boundary condition then becomes 



Z PP 



n cosh (pjt) 



8=1 



Q Q 

fii,i + fii, + ^o,i - s S n ( —q^ ) ^om 



Furthermore, the free energy density per ksT of the system defined by 



then takes the form 

f pp = - ^ N \~ ^ hi 2 - ]T In [cosh {(3 J;)] 



i=l 



(69) 



(70) 



(71) 



Periodic- antiperiodic boundary condition 
For pa boundary condition, equation ( p8|) is replaced by 
B" (CM) + = - [(CM) + + (CA) + B~ - (CAY B + + {CAY B~] , 
and the partition function has the form 

-1 N. 



r,pa _ _2 N 3 

" 2 



rib 



J^J cosh (pji 



i=l 



S7i 1 + Jl 



i,0 "0, 5 



fi n 1 + sgn ( ) fi o 



0, 



(72) 



(73) 



where the superscript a refers to antiperiodic boundary condition. The corresponding 
free energy density per k B T is 

(N s - 1) 



/ 



pa 



N s 



ln2-^m [cosh (f3Ji)] 



i=l 



Q Q 

fii i + fii - fi i + sgn ( — — - ) fi o 



(74) 



3.3. Antiperiodic-periodic boundary condition 

Similarly, for ap boundary condition, equation fl3"8| ) is replaced by 

B + {CAY = - [{CAY B + - (CA) + B~ + {CAY B+ + (CAY B~] 



and 



7ap 



Y[ cosh {(3Ji 



i=l 



fii 1 — fii Q + fi 1 + sgn 



6-6 c 
~Q7 



fi 



00 



(75) 



(76) 
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/ 



i=l 



fiii - % o + ^o,i + s S n ( ) fi oo 



(77) 



3.4- Antiperiodic-antiperiodic boundary condition 



For aa boundary condition, equation ( |38"D becomes 

B~ (CA)~ = - [- (CA) + B + + (CA) + B- + (CA)' B + + (CA)~ B~] 

and the partition function is 



Z aa _ _ 2 

2 



1=1 



-Qi i +Qi 



fi i + sgn 



e-e c 



n 



00 



(78) 



(79) 



The corresponding free energy density per ksT is 

(N s - 1) 

m z — 

i=l 



r 



ln2-^Sn [cosh (/3J;)] 



(80) 



1 1/2 



1) 



Note that by taking t 3 = 0, n b = 2 and N s = N b = MN, we have A 3 = 

AT-l iV-1 

^=nil [a -A iC os 2 -^-A 2 cos 2 -^ 

and all the results we obtained reduce to those of sq lattice. 

Accordingly, the critical temperature can be determined in the thermodynamic 
limit from the zero of the free energy contributed by the zero mode, 

A - A 1 - A 2 - A 3 = 0. (82) 

It follows that for isotropic coupling, we have 



0, 



im(i + V5 



-1 



2.269185... 



for sq lattice with 9 = ksT / J, 



0c = 

for pt lattice and 

9c = 

for he lattice. 



Im(V3 



3.640956. 



iln( 2 + V3 



-1 



1.518651. 



(83) 



55) 
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4. Specific heat 

The specific heat per spin C/ks for the Ising model onMxAf sq, pt and he lattices with 
isotropic couplings are shown, respectively, in figure 2(a), 3(a), 4(a) for M/N = 1, and 
in figure 2(b), 3(b), 4(b) for M/N = 1/2. Figure 3(c) and 4(c) show, respectively, results 
for pt and he lattices under pa and aa boundary conditions and for M/N = 1, 1/2, 1/4. 
In general, for three lattices with the same lattice size, the specific heat under pp 
boundary condition is always larger than those under other boundary conditions. Note 
that for sq lattices with M/N = 1, C pa and C aa are distinct in figure 2(a), but for 
pt and he lattices with M/N = 1 in figure 3(a) and 4(a), they coincide and are non- 
distinguishable due to the last term in the bracket of equation (|68"D , which is associated 
with the structure symmetries of pt and he lattices. These behaviours can be violated 
by taking aspect ratio £ = M/N ^ 1, and the results are shown in figures 3(b), (c) and 
4(b), (c). 

We further study the displacements of the maximum of C pp and C pa . The shift 
behaviours of the maximum in Cjvjv (T) are shown in figure 5. The slopes of the curves 
implies the rates of approach of C pp and C pa to their limiting behaviours. For periodic- 
periodic boundary condition, these lattices have linear behaviours in iV — > oo and can 



be described by the formula |36|| 



(T c ^max) O A7 - (or\ 

— ~ — , as iV — > oo. 86 

T c N y ' 

For periodic-antiperiodic boundary condition, the corresponding formula is also provided 
by finite-size scaling ansatz. However, for numerical analysis, instead of equation (|86|), 
we use 

(T c — T max ) a bi 62 

As a result, we have a pp = 0.360, 6ff s = -0.47, a pa = 0.18, b\ a s = -2.19, for the sq 
lattice, a pp = 0.363, = -0.91, of = 0.09, b^ t = 0.60 for the pt lattice, a p h p = 0.268, 
h V \h = 0.24, a^ a = 0.09, = 0.87 for the he lattice, and the value of b 2 is of the 
order of 1. The values of a pp is larger than a pa for three lattices, and this implies the 
approach to limiting behaviour for pp boundary condition is faster than pa boundary 
conditions. Since the logarithmic divergence of the specific heat is independent of 
boundary conditions and can not be used to distinguish C pp and C pa of large lattice, 
then the values of a may be used to distinguish two boundary conditions. 



5. Discussion 



We have solved the exact partition functions of M x pt and he lattices with different 
boundary conditions. These results can provide the analytical background for further 
studies on the effects of lattice structures and boundary conditions on critical properties 
and critical finite-size corrections of the Ising model. 
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Firstly, universal finite-size scaling functions for critical systems have received much 
attention in recent years |L5], |16|, [T7|, |2T], [J7|, [38| , and it is well known that the finite- 
size scaling functions depend on the boundary conditions \T%\. Hu, Lin and Chen, and 
Okabe and Kikuchi have discussed the difference in the finite-size scaling functions for 
lattice models under periodic boundary and free boundary conditions in connection with 
the universal finite-size scaling function for percolation problem and Ising model (TJJ . 
respectively. Other boundary conditions, such as the Ising model on an M x N simple- 
quartic lattice embedded on a Mobius strip and Klein bottle also has been studied PH . 
Kaneda and Okabe found that there is interesting aspect ratio dependence of the value 
of the Binder parameter at criticality for various boundary conditions |19| . It is then 
interesting to have a rigorous test of finite-size scaling function and critical finite-size 
corrections for different planar Ising model under various boundaries. 

In addition, by using Monte Carlo method, Hu, Lin and Chen [15, [16|], Tomita, 
Okabe and Hu |2T| have found that the universal finite-size scaling functions of the 
scaled quantities for sq, pt and he lattices depend on the aspect ratios and have very 
good universal finite-size scaling behaviours when the aspect ratios of these lattices have 
the proportions 1 : y/3/2 : a/3. This further implies lattice-structure-dependence of the 
universal finite-size scaling function and it would be a rigorous test from analytical 
aspect. To aforementioned topics, we have found finite-size scaling behaviours for sq, pt 
and he lattices under period-aperiodic boundary conditions. By selecting a very small 
numbers of nonuniversal metric factors, we have further found very good universal finite- 
size scaling behaviours for these lattices, and the results will be presented in other paper. 

Finally, the discussion of the specific heat in this paper also inspire another problem. 
Quite recently, Izmailian and Hu have found exact amplitude ratio and finite-size 
corrections for the M x N sq lattice Ising mode on a torus M, and new sets of the 
universal amplitude ratios of subdominant correction to scaling amplitudes ||. The 
results of section 4 suggest that a pp /a pa for sq, pt and he lattices are roughly 2, 4, and 
3. The question is "is there exact relations between a pp and a pa for these lattices ?" It 
is interesting to study this question and to have a heuristic argument for this simple 
relation. 
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Figure 1. (a)The global structure of the triangular lattice used in this paper. A basic 
cell of the lattice site is given by (m, n), and the coupling constants are Ji, J2 and J3. 
(b)The global structure of the honeycomb lattice used in this paper. Each basic cell 
contains an inner Ising spin uq. 



Figure 2. The specific heat per spin for (a) N x N square Ising lattices with isotropic 
couplings under pp, pa, ap and aa boundary conditions, and (b) M x N square Ising 
lattices with isotropic couplings and aspect ratio M/N — 1/2 under pp, pa, ap and aa 
boundary conditions. The critical point 6 C is marked by a vertical line. 



Figure 3. The specific heat per spin for (a) N x N plane-triangular Ising lattices 
with isotropic couplings under pp, pa, ap and aa boundary conditions, (b) M x N 
plane-triangular Ising lattices with isotropic couplings and aspect ratio M/N =1/2 
under pp, pa, ap and aa boundary conditions, and (c) M x N plane-triangular Ising 
lattices with isotropic couplings and aspect ratio M/N = 1, 1/2, 1/4 under pa and aa 
boundary conditions. The critical point 6 C is marked by a vertical line. 



Figure 4. The specific heat per spin for (a) N x TV honeycomb Ising lattices with 
isotropic couplings under pp, pa, ap and aa boundary conditions, (b) MxN honeycomb 
Ising lattices with isotropic couplings and aspect ratio M/N — 1/2 under pp, pa, ap 
and aa boundary conditions, and (c) M x JV honeycomb Ising lattices with isotropic 
couplings and aspect ratio M/N = 1, 1/2, 1/4 under pa and aa boundary conditions. 
The critical point 9 C is marked by a vertical line. 



Figure 5. (a) Variation of (Tmax — T c ) with finite N for N x N square Ising lattices 
with isotropic couplings under pp and pa boundary conditions. The broken lines are 
given by (T max — T c ) /T c — a/N and indicate the limiting behaviour as N — > 00. (b) 
Variation of (T max — T c ) with finite N for N x N plane-triangular Ising lattices with 
isotropic couplings under pp and pa boundary conditions, (c) Variation of (T max — T c ) 
with finite N for N x N honeycomb Ising lattices with isotropic couplings under pp 
and pa boundary conditions. 



